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ABSTRACT 

We study an acceleration of electrons and positrons in the relativistic magnetic 
field reconnection using a 2.5-D particle-in-cell electromagnetic relativistic code. 
We consider the model with two current sheets and periodic boundary conditions. 
The electrons and positrons are very effectively accelerated during the tearing and 
coalescence processes of the reconnection. We found that near the X-points of 
the reconnection the positions of electrons and positrons differ. This separation 
process is in agreement with those studied in the previous papers analytically 
or by test particle simulations. We expect that in dependence on the magnetic 
field connectivity this local separation can lead to global spatial separation of the 
accelerated electrons and positrons. A similar simulation in the electron-proton 
plasma with the proton-electron mass ratio mi/m e = 16 is made. 

Subject headings: Acceleration of particles - Plasmas - Relativity 



1. INTRODUCTION 



Magnetic reconnection is the key process in conversion of the magnetic field energy 
into particle kinetic energy. It is well accepted that it plays a crucial role in the Earth's 
magnetotail, solar flares and accretion discs (Priest & Forbes 2000, Drake et al. 2005, 
Pritchett 2006). The relativistic reconnection in electron-positron plasmas is proposed for 
high energy astrophysical phenomena including the jets from active galactic nuclei (Lesch & 
Birk 1998, Larrabee et al. 2003, Wardle et al. 1998), pulsar winds (Coroniti 1990, Michel 
1994, Lyubarsky & Kirk 2001) and models of gamma-ray bursts (Drenkhahn 2002a,b). The 
relativistic reconnection and particle acceleration in pair plasmas was studied numerically for 
the first time by Zenitani & Hoshino (2001, 2005). The effectiveness of such an acceleration 
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and corresponding synchrotron spectra has been computed in detail in the paper by Jaroschek 
et al. (2004a,b). Bessho & Bhattacharjee (2005) have shown that this fast reconnection is 
caused by the off-diagonal components of the pressure tensor. 

Recently, the RHESS1 observations of the solar flare on 2002 July 23 have revealed a 
separation of the gamma-ray source from any of those observed in the hard X-ray emission. 
It has been interpreted as a spatial separation of energetic electrons and protons (Share et al. 
2003). Based on the analytical and test particle approach Zharkova & Gordovskyy (2004) 
have explained this separation by the asymmetry in acceleration of electrons and protons in 
the reconnecting non-neutral current sheet, see also the papers by Martens & Young (1990), 
Zhu & Parks (1993), Litvinenko (1996). 

This separation acceleration can be even more distinct in pair plasmas due to the same 
mass of electrons and positrons. Therefore in this paper using the particle-in-cell modelling 
we study this process in the electron-positron plasma in detail. 

2. MODEL 

We used a 2.5-D (2D3V - 2 spatial and 3 velocity components) fully relativistic electro- 
magnetic particle-in-cell code (Saito & Sakai 2004). The system size is L x x L y = 2000 A x 
600A = 200<i e x 60d e , where A (=1) is a grid size, d e = c/u pe is the electron inertial length, 
c is the speed of light and u pe is the plasma frequency. 

Two 2-D current sheets with the guiding magnetic field B 2 are initiated along the lines 
y = 150A and y = 450A. The periodic boundary conditions are used. The half-width of 
both the current sheets is 10A = d e . The initial magnetic field is (see also Karlicky & Barta 
2007) 

B = (B x , By, B z ), 
B x = -B for y < 140A, 
B x = (y- 150)£ /10 for 140A < y < 160A, 
B x = B for 160A < y < 440A, 
B x — —(y — 450) S / 10 for 440 A < y< 460A, 

B x = -B for y > 460A, 
By — 0, B z — Bq. 

We consider the electron-positron plasma. In each numerical cell located out of the 
current sheet we initiated n = 60 electrons and n = 60 positrons. In this region out of 
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the current sheet we define the plasma frequency for the time unit oC, 1 . The time step in 
computations is co> pe At = 0.05. The total amount of particles in the model is 172 millions. The 
initial number density is enhanced in the current sheets just to keep the pressure equilibrium 
in the current sheet. The particle distribution is taken as by Zenitani & Hoshino (2001) in the 
form of f ~ exp —m[u 2 x + u 2 + {u z — U) 2 ]/2T, where the velocity u is related to the particle 
velocity v as u = v 7 = v/[l — (v / c) 2 ] 1 / 2 , U is the drift velocity, m is the electron rest mass, 
T is the plasma temperature and c is the speed of light. The mean initial thermal energy 
of electrons and positrons is taken the same as 0.45 mc 2 . We neglect any collisions, pair 
production, and pair annihilation of pair plasmas. 

Due to our interest about the reconnection processes in the relativistic plasma with 
high magnetic field we consider cases with low-/? plasmas. The plasma beta parameter and 
the ratio of the electron-cyclotron and electron-plasma frequencies in the region out of the 
current sheets are chosen as /?=0.11, u> ce /u} pe = 4 (Case I) and (3=0.5, u ce /u pe = 1.9 (Case 
II). For comparison one run was made for the parameters as in Case I, but without the 
guiding magnetic field, i.e. B z = 0. 

Furthermore the same processes are modelled in the electron-proton plasma (Case III) 
with the proton-electron mass ratio mi/m e = 16. The proton and electron temperature is 
taken the same Tj = T e . The parameters are (3 = 0.11 and uj ce /uj pe = 4. The mean initial 
thermal energy is 0.45 mc 2 . 

All computations were performed on the parallel computer OCAS (Ondfejov Cluster 



for Astrophysical Simulations), see http://wave.asu.cas.cz/ocas 



3. RESULTS 

Due to the tearing mode instability the current sheet tears into O-type islands (plas- 
moids) which later on coalesce into larger ones. During these processes both the electrons 
and positrons are accelerated. Figure 1 shows an evolution of the electron distribution func- 
tion f(E) (for Case I) in dependence on the energy E in the whole numerical plane (upper 
part of Figure 1) as well as in one selected location (bottom part of Figure 1, for the location 
see Figure 4), where much harder spectrum can be seen. While the spectral index at u pe t = 
600 in the whole plane is -3.3, in the selected location is about -1.9. There are other such 
places, especially at the locations where high-energy electrons and positrons are produced, 
see Figure 4. This acceleration process is very efficient and looks to have a 2-step character 
in time as shown in Figure 2, where an evolution of the number of accelerated electrons in 
the ratio to the total number of electrons (in percents) with the energy E/m 2 = 7 > 4 is 
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shown. An analysis shows that the first step of the acceleration (up to about u pe t = 400) is 
connected with tearing processes and in the second step (after u pe t = 400) the main process 
is the coalescence of plasmoids. The acceleration with the higher value of u ce /uj pe = 4 (full 
line) (Case I) is more efficient than that with uj ce /uj pe = 1.9 (dashed line) (Case II) (Figure 
2). 

Similarly as in previous studies (Drake et al. 2005, Pritchett 2006, Jaroschek et al. 
2004a,b) the electrons and positrons are accelerated in the electric field near the X-points 
formed during the tearing and coalescence processes. But we found that during this ac- 
celeration process (Figures 3 and 4) the electrons and positrons are moving into different 
locations around the O-type magnetic structures (plasmoids) and thus they are spatially 
separated (e.g. see the region around x = 500 A and y = 450 A in Figure 4). To understand 
this separation process we analyzed the electric field near the X-point of the reconnection 
(Figures 5 and 6). As seen here the electrons (the asterisks) are located at the borders of 
the areas with the enhanced (—E y and — E x ) and (+E y and +E X ). On the other hand, the 
positrons are located along the remaining two borders (see Figure 3, bottom part). The 
electric component E z along the line y = 150A is negative between the O-type magnetic 
structures (Figure 6). But only near the X-point of the reconnection this electric field (E z ) 
deviates from that of the inductive one — v x B/c (where v is the plasma velocity), see Fig- 
ure 7, where their profiles are shown at two times {oj pe t = 50 and 100). This deviation defines 
the diffusion region of the reconnection. As concerns the magnetic field in the early stage 
of the reconnection, the magnetic field component B y is positive in the region x > 1000A 
and negative for x < 1000A. The structure of the magnetic field together with the plasma 
velocity pattern in the regions from the X-point to the magnetic island centers resembles 
to that of the collapsing magnetic trap (Giuliani et al. 2005; Karlicky & Barta 2006). In 
such a structure the particles are also accelerated, but not separated as found here near the 
X-points. This additional acceleration process is known as that in the contracting magnetic 
islands (Drake et al. 2006). 

Furthermore, we made the similar computations, but without the guiding magnetic field 
component (B z = 0). In this case no separation of electrons and positrons was found. Also 
the electric field structure was different from that presented in Figures 5 and 6. 

Finally, we made similar computations also for the electron-" proton" plasma, with the 
electron-proton ratio mi/m e = 16 (Case III). Similarly as in the previous cases, Figure 8 
shows that accelerated electrons and protons move to different positions. Comparing this 
Case III (dotted line in Figure 2) with Case I (full line) the number of accelerated electrons 
is reduced. 
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4. DISCUSSION AND CONCLUSIONS 

The present simulations show that the magnetic reconnection with the guiding mag- 
netic field accelerate electrons to different positions around the plasmoid than for positrons 
or "protons". If the magnetic field connectivity (in the z-direction) from the upper and 
bottom part of the plasmoid differs then the accelerated electrons and positrons (or protons) 
move into quite different locations as observed by RHESSI. The separation of particles with 
different electric charge is a natural consequence of the acceleration in direct electric field 
near the X-points of the magnetic field reconnection. 

In agreement with Litvinenko (1996) and Zharkova & Gordovskyy (2004) we found that 
the separation process is due to a presence of the non-zero guiding magnetic field (non- 
neutral current sheet). Namely, our simulations with the zero guiding magnetic field show 
no such separations. 

Zharkova & Gordovskyy (2004) have shown that the separation direction (in the present 
designation of the electric and magnetic field components) depends on the sign of the term 
q 3 B y B z E z , where q is the electron (-e) or positron (+e) charge (see the relation (8) in their 
paper). Considering the direction of the magnetic and electric fields in our case (Figures 5 
and 6) it can be shown the separation direction found agrees to this relation. This result 
agrees also to the relations presented in the paper by Litvinenko (1996). 

The similar separation process is found also for the reconnection in the "proton" -electron 
plasma with the proton-electron mass ratio rrii/m e = 16. This mass ratio is not realistic and 
is taken due to computer limitations. Nevertheless, we expect that such a separation process 
will be confirmed by future computations also for the real proton-electron mass ratio. 

Comparing the acceleration process for the electron-positron (full line) and electron- 
proton (mi/m e — 16, dotted line) plasma in Figure 2, we found that N m ./ me=1 /N m ./ me=16 is 
1.18, i.e. the number of accelerated electrons N depends on the proton-electron mass ratio 
as Ni/N m ./ me ~ l/(mj/m e ) ' 0625 . If this relation is valid also for the real proton-electron 
mass ratio then we can write Ni/Ni 8 38 = 1-6, which gives enough accelerated electrons also 
for the real electron-proton plasma. 

Comparing the present modelling with previous studies the most similar simulation is 
that of Zenitani & Hoshino (2001), especially due to initial high thermal plasma energy. But 
in their model no guiding magnetic field, which is crucial for the particle separation, is con- 
sidered. Furthermore, contrary to our start from noise level they initiate the reconnection 
by magnetic field perturbation which can influence the separation process, too. The maxi- 
mum energies of accelerated electrons in both models are comparable. But in our model the 
reconnection process is about three times slower than that in Zenitani & Hoshino (2001). 
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We think that it is due to the initial magnetic field perturbation. 

Although the acceleration in the contracting magnetic islands does not separate particles 
of opposite electric charges this process is important for global acceleration. But this process 
is even more complicated than presented in the paper by Drake et al. (2006). Namely, not 
only parallel energy of particles increases due to reflection from the ends of contracting 
magnetic islands (as described by the relation (1) in Drake et al. (2006)) but also the 
perpendicular energy of particles E± can increase due to the betatron type of the acceleration, 
which follows from the conservation of the magnetic moment fi = E±/B in the region 
with the increasing magnetic field B (see Karlicky & Barta 2006). In our simulations the 
contracting acceleration is time varying, therefore let us compare its efficiency with that of 
the acceleration near the X-point at one specific time. Using the relation (1) of Drake et al. 
(2006) we derived the electric field equivalent to this process as E eq = (v x B 2 , v m e )/(8 x e 
B 2 ), where 25 x is the length of the island, B x and B are the reconnecting and total magnetic 
fields, v x is the contracting velocity, v is the mean electron velocity, and m e is the electron 
mass. For the parameters in Case I at u pe t = 100 it gives the equivalent electric field in the 
contracting magnetic island E eq one order of magnitude lower than the electric field E z at the 
X-point region. It means that at this moment the acceleration near the X-point dominates 
over that in the contracting magnetic islands. 

This research was supported by the Centre for Theoretical Astrophysics and by Grants 
IAA300030701 of the Academy of Sciences of the Czech Republic. M. K. thanks Prof. J. I. 
Sakai and Dr. S. Saito for many useful discussions concerning numerical modelling. Author 
also thanks to the referee for constructive comments improving the paper. 
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Fig. I. — The energy distribution of electrons at the initial state (dotted line), at tu pe t = 600 
(dashed line), and at uo pe t = 1000 (full line) (Case I). The upper part: The distributions in 
the whole numerical plane. The bottom part: The distributions in the selected location, in 
the circle centered at x=550A and y=450A with the radius r = 50A. For this location, see 
Figure 4. 




Fig. 2. — The number of accelerated electrons with E/mc 2 > 4 expressed in the ratio to 
total electron number for Case I (full line), Case II (dashed line), and Case III (dotted line). 
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Fig. 3. — The spatial distribution of accelerated electrons (dots in upper part) and positrons 
(dots in bottom part) with the energy E/mc 2 > 10 superposed on the magnetic field lines 
projected to the x — y plane at u pe t = 150 (Case I). Compare locations of electrons and 
positrons. 




Fig. 4. — The spatial distribution of accelerated electrons (dots in the upper part) and 
positrons (dots in the bottom part) with the energy E/mc 2 > 13 superposed on the magnetic 
field lines projected to the x — y plane at u pe t = 600 (Case I). For the separation of electrons 
and positrons, see the region around x = 500A and y = 450A. 
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Fig. 5. — The spatial distributions of the electric field components (E/Bo) around the X- 
point in the current sheet near y = 150A at uo pe t = 100 (Case I), compare with Figure 6. 
The asterisks mean accelerated electrons for the energy E/mc 2 > 10. The thick full contour 
means the area of enhanced | E z | component (level E z j Bo— -0.22). (Remark: In the upper 
current sheet with oppositely oriented current along the line y=450A the E z component is 
oppositely oriented, i.e. positive.) The thick dashed contour is the area of the enhanced 
E y (level E y /B — 0.36). The thick dotted contour means the area with oppositely oriented 
E y (level E y f Bo— -0.36). The thin dashed contour is the area of the enhanced E x (level 
E x /B o —0. 07). The thin dotted contour means the area with oppositely oriented E x (level 
E x /B =-0.07). 




Fig. 6. — The electric field components (E x (dashed line), E y (dotted line), and E z (full 
line)) along the center of the current sheet, i.e. along the line |/=150A at u pe t = 100 (Case 
I) (compare with Figure 5). 




Fig. 7. — The electric field component E z (full line) and the term — v x B/c (dashed line) 
at two times (Case I): at u pe t = 50 (upper part), and at ui pe t = 100 (bottom part). 
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Fig. 8. — The spatial distribution of accelerated electrons (dots in the upper part) with 
the energy E/m e (? > 5 and "protons" (dots in the bottom part) with E/rriiC 2 > 1.4 near 
the X-point of the magnetic field structure in the x — y plane at u pe t = 50 (Case III). The 
proton-electron mass ratio is 16. 



